Linear magnetoresistance from Dirac-like fermions in graphite 
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We show that magnetoresistance of Bernal-stacked graphite (with the magnetic field B parallel to 
the c-axis and the current in the ab plane ) scales linearly with the magnetic field over an interval of 
classically weak fields. The linearity is related to the presence of extremely light, Dirac-like carriers 
near the H (H')~ points of the Brillouin zone. The Hall resistivity in this interval also shows a 
non-analytic, B\n\B\ behavior, and is dominated by holes. 
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I. INTRODUCTION 

Dirac fermions have become a focus of great interest 
in condensed matter physics in recent years. Being mass- 
less, Dirac fermions respond to the magnetic field much 
more strongly than their massive counterparts. A com- 
mon platform for studying Dirac fermions is graphenc, 
where their sensitivity to the magnetic field allows for 
the observation of the quantum Hall effect even at room 
temperature^ Stacking up graphene layers in a Bernal 
(ABAB...) way, one forms a three-dimensional crystal 
of graphite. For the most part of the Brillouin zone 
(BZ) and at energies smaller than interlayer hopping 
(~ 0.4 eV), the charge carriers in graphite have little 
resemblance to their Dirac "ancestors": they are just 
massive electrons and holes. Graphenic ancestry, how- 
ever, makes the effective mass depends on the momen- 
tum along the c axis in such a way that it vanishes at the 
top and bottom edges of the BZ (H(H')- points). Tiny 
regions around these points harbor Dirac-like fermions 
similar to those in a single-layer graphene. Landau levels 
of these Dirac fermions in graphite have been observed 
by tunneling^ and opticaP spectroscopies. Their contri- 
bution to transport, however, is not clear; in particular, 
their role in quantum magnetoscillations has been a su b- 
ject of recent discussion and remains controversialpEl 

Although transport in graphite has been well studied 
in the past, there remain quite a few open questions. For 
example, neither the magnitude nor the temperature or 
magnetic- field dependences of the c-axis conductivity can 
be satisfactorily explained within the Boltzmann trans- 
port theoryPtlSl i n comparison, the in-plane conductivity 
is relatively better understood; for example, its tempera- 
ture dependence is well explained using the scmiclassical 
Boltzmann equation! 11 " 14 ! The success, however, is only 
partial as the explanation of field-dependence of the in- 
plane magnetoresistance presents some difficulty. Experi- 
mentally, the in-plane magnetoresistance (MR) (with the 
magnetic field along the c-axis and the electrical current 
in the afc-plane) is often found to depend linearly on the 
magnetic fielcpSHm which, at first glance, seems to con- 
tradict transport theory. Even more surprisingly, linear 
MR spans over a wide range of fields, beginning at classi- 
cally weak fields and persisting up to the ultra-quantum 



regime and beyondP Although the linearity in the ultra- 
quantum regime can be accounted for by taking into ac- 
count the field dependence of the scattering timeP the 
linear behavior in semi-classical fields still lacks a proper 
understanding. Although it is sometimes ascribed to ex- 
trinsic reasons, such as macroscopic inhomogeneities, the 
issue is far from being settled.^ 

Since the magnetic field is a (pseudo) vector, MR 
can only be a function of B 2 . In general, MR behaves 
quadratically for weak fields and either saturates, if elec- 
tron orbits are closed, or grows quadratically in strong 
fieldsp^ED if orbits are open or in compensated metals. A 
linear and thus non-analytic dependence on the field indi- 
cates some non-trivial physics. Experimentally, graphite 
is not the only material that exhibits linear MR — it is 
found to occur in many other materials as wellpSl yet de- 
tailed understanding of this effect exists only in a hand- 
ful of situations.^ For example, it has been shown that a 
non-analytic dependence on the fi eld ca n arise due to spe- 
cial features of the Fermi surface J22123I w hich makes one 
wonder whether graphite has any such peculiarities too. 
In fact, previous studies of MR within the Boltzmann 
transport theory have already hinted at such a possibil- 
ity: numerical calculations have shown that at low tem- 
peratures and in weak fields, MR behaves as B n with 
n < 2, which indicates a departure from the canonical 
behavior. 

In this paper, we show, using a simplified yet con- 
sistent model for the energy spectrum of charge carri- 
ers and semiclassical Boltzmann equation, that linear 
MR is an inherent property of graphite bandstructure. 
The linearity stems from very light, Dirac-like carri- 
ers near the H(H')- points of the BZ. The origin of 
the effect can be understood without detailed calcula- 
tions. Indeed, anywhere but in the immediate vicini- 
ties of the H(H')- points, electron-and hole-like car- 
riers in graphite move in the a6-plane as free parti- 
cles but with a mass that depends on the momentum 
along the c-axis (k z ). The corresponding spectrum is 
described by £k = /(&«) ± k 2 /2m*(k z ), where k p is 
the in-plane momentum, f(k z ) is some function of k z 
and m*(k z ) ex cos(k z c/2) with c being the c-axis lat- 
tice constant. The in-plane magnetoconductivity con- 
tains the usual factor [l + uj 2 (k z )T 2 ] averaged over k z 



(here ui c (k z ) = eB/m*(k z ) and r is the relaxation time). 
Near the H(H')- points (k z — ±n/c), the effective mass 
vanishes as n/c — k z , the cyclotron frequency diverges, 
and the integral over k z behaves as \B\ instead of B 2 . 
Likewise, the Hall conductivity contains an average of 
u c {k z ) [l + ui 2 (k z )T 2 ] , which also diverges, albeit only 
logarithmically, near the H(H')- points. As a result, 
the Hall conductivity behaves as i?ln|_B|. The width 
of the regions contributing to non-analytic MR shrinks 
with B. For B below some -Bo, which is still parametri- 
cally smaller than the weak-to-strong field crossover B c 
(defined such that eB c T/m*(0) = 1), these regions be- 
come so narrow that the quadratic model of the spectrum 
ceases to be valid. Using the full spectrum instead of the 
quadratic approximation, one can show that MR is ana- 
lytic for B <C B a . In classically strong fields (B ^> B c ), 
MR is again analytic. Thus, non-analytic MR occurs in 
intermediate fields, i.e., for Bo C B <C B c . An order 
of magnitude estimate for B is (|£f|/7i)-B c ps 0.06B C , 
where Ef is the Fermi energy and 71 is the nearest-plane 
hopping; thus the interval of fields between B c and Bq 
is sufficiently wide. In experiment, however, linear MR 
often spans the entire interval of magnetic field: from the 
weak-field regime to the ultra-quantum limit. Therefore, 
intrinsic linear MR, predicted in this paper, cannot ex- 
plain the data for all magnetic fields. It is possible that 
other factors, such as macroscopic inhomogeneities, are 
responsible for linear MR in strong magnetic field. 

The rest of the paper is organized as follows. Sec.[lT|de- 
scribes the model for the electronic spectrum of graphite 
which allows for an analytic calculation of the conduc- 
tivity. In Sec. |III[ we calculate the field dependences 
of the magnetoconductivity, first in the quadratic- "non- 
relativistic" approximation (Sec. Ill A I, and then for the 
full spectrum (Sec. Ill B ). The effect of macroscopic inho- 
mogeneities is discussed in Sec. HIC 1 An issue whether 



Dirac-like carriers in graphite play a role in quantum os- 
cillations is addressed in Sec. lHIC"2l Our conclusions are 
given in Sec. |IV| 
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FIG. 1. Schematic representation of the Fermi surface of 
graphite. Non-relativistic massive fermions occupy most of 
the Fermi surface, except for the narrow regions near the 
H(H')- points, occupied by massless Dirac fermions. Here 
#0 = cos _1 \/ej?/272 ~ 0.62 marks the point at which the 
Fermi energy intersects the electron and hole bands. 



ergy spectra of the conduction and the valence bands, 
respectively, can be written as 
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II. 



ENERGY SPECTRUM OF GRAPHITE 



Ideal graphite consists of graphene planes stacked on 
top of each other in the Bernal way (ABAB...). Its band 
structure is usually described by the Slonczewski-Weiss- 
McClure (SWMc) modeP characterized by seven param- 
eters: 70 ■ • ■ 75 and A. Here, 70 and 71 denote the in- 
plane and out-of-plane nearest-neighbor hopping terms, 
correspondingly; 72 ■ • ■ 75 denote various next-nearest- 
neighbor hopping terms; and A arises due to the differ- 
ence between the on-site energies of the A and B carbon 
atoms. In terms of energy scales P 70 ps 3.2 eV is the 
largest one, followed by 71 ps 0.4 eV, 73 ps 0.3 eV, and 
74 w 0.1 eV, followed by 72 « -0.02 eV, 75 ps 0.01 eV, 
and A ps 0.01 eV. A closed form of the energy spectrum 
can be obtained only if 73 is put to zero. Under this 
approximation (discussed in more detail below), the en- 



where e° 1>2 = A± 7l r+ i 75 r 2 , e° 3 = i 72 r 2 , v p = ^7o«, 
and a(c) is the in- (out-of-) plane lattice constant. 

To facilitate an analytical calculation of the conductiv- 
ity later, we simplify these expressions as much as pos- 
sible, while keeping the physical content intact. First, 
we note that the parameter 74, which enters as a ratio 
74/70 ~ 0.1, only introduces a weak k z dispersion in the 
in-plane velocity and thus can be neglected. Next, al- 
though the parameters 72, 75, and 76 are of the same 
order of magnitude, their roles in the spectrum are very 
different. Namely, 72 plays a crucial role as it determines 
the band overlap, and hence small («3x 10 18 cm~ 3 as 
T — > 0) but non-zero carrier concentration. On the other 
hand, 75 and A do not lead to qualitative changes in the 
spectrum and thus can be neglected. With these simpli- 
fications, the minimal model for the energy spectrum in 



graphite can be written as 
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(For 72 = 0, Eq. ([2| reduces to the original Wallace 
modePS which describes graphite as a zero-gap semicon- 
ductor.) We will refer to the spectrum in Eq. pi) as to 
"relativistic" . Equation pi) describes two groups of mas- 
sive Dirac fermions with the "rest mass" which vanishes 
at the H(H')- points k z = ±7r/c, where T = (cf. Fig. 
[TJ). Away from the H{H')- points one can expand Eq. (fTl) 
in fc p , which gives a "non- relativistic" spectrum with the 
effective mass varying with k z : 



where (f> is the azimuthal angle, and all the smaller band- 
structure parameters were neglected in the same way 
as in Eq. ([3]). Not too close to the H{H')- points, 
i.e., for |r| ~ 1, the trigonal corrections are smaller, 
though not in order of magnitude, than the first term 
(2) [kp/2m*(k z ) ~ sp]. At the K point, for example, the 
second term amounts to ~ 73/7ov7i/£f£f ~ 0.45e^ 

and the third term to ~ (73/70 \/ii/£f) 2 £f ~ 0.23£f, 
so that the inclusion of 73 does lead to non-negligible 
changes in the spectrum. However, since the trigonal cor- 
rections are proportional to a power of T, they become 
much smaller near the H(H')- points. Therefore, in the 
regime when the magnetoconductivity is dominated by 
the carriers near the H(H') points, 73 can indeed be ne- 
glected. This is the regime we will focus on in the next 
Section. 
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The Fermi energy, ep , is determined by a balance of the 
two terms in Eq. pi , each of which is of order ep . There- 
fore, we did not neglect 72 in the first term, as this would 
have rendered ep to zero, but neglected it in the second 
term, where it it would have only given a small correction 
to the effective mass. 

Following the tradition, we will refer to the ± branches 
of the spectrum as electron/hole bands, although in fact 
the bands have a mixed electron-hole character, depend- 
ing on the magntitude and direction of k. For example, 
the effective mass for in-plane motion is positive (nega- 
tive) for the +(— ) band, hence the + (— ) band corre- 
sponds to in- plane electrons (holes). On the other hand, 
the two bands are degenerate at k p = 0, while the effec- 
tive mass along the c axis is positive for \k z \ < tt/2c and 
negative for 7r/2c < \k z \ < tt/c (recall that 72 < 0). 

In the non-relativistic approximation, the electroneu- 
trality condition fixes the Fermi energy to Ep = 
(4/3)72 ~ —26 meV, very close to the commonly ac- 
cepted value of —25 meVP Massless Dirac fermions with 
dispersions £k = ~v p k p , located near the H(H')- points, 
do not modify significantly the charge balance. However, 
as we will show in the next Section, the vanishing mass of 
the in-plane motion near the H(H')- points affects dra- 
matically both the longitudinal and Hall conductivities. 

Before moving on to the calculation of the conductiv- 
ities, we need to justify the neglecting 73 in Eq. ([I]). A 
perturbat ion theory in 73 results in a trigonally warped 
spectrumP25l 
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III. 



LINEAR MAGNETORESISTANCE 



A. Magnetoconductivity for the non-relativistic 
energy spectrum 

We now calculate the components of the magnetocon- 
ductivity tensor in the non-relativistic approximation for 
the energy spectrum described by Eq. ([3|. We use the 
linearized Boltzmann equation in the relaxation time ap- 
proximation, 
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df° 
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de 



e(v x B) 
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5(k), (5) 



where e > is the magnitude of the electron charge, f° 
is the Fermi function, and g(k) is the non-equilibrium 
part of the distribution function. In general, the Boltz- 
mann equation cannot be solved in a closed form for 
an arbitrary spectrum (even within the relaxation-time 
approximation). However, since Eq. ([3]) is isotropic in 
the in-plane direction, one can solve Eq. (pi) exactly for 
the case of the magnetic field along the normal to the 
plane. Assume that g(k) = v ■ A, where A is an in- 
plane vector represented in terms of E and B x E as 
A = aE + 6B x EP^l Then g(k) contains only the in- 
plane velocity v p = k p /?7i* (k z ). A simplifying feature of 
the problem with in-plane isotropy is that the coefficients 
a and b are allowed to be functions only of k z but not of 
k p . Substituting <?(k) into Eq. (15J) and solving for a and 
b, we obtain 



A = re 
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The diagonal components of the magnetoconductivity for 
each of the electron and hole bands are given by 



4e 2 r f vl ( df° 



ot(B) = <4(B) = ^3 J 1 + ^ {kz)T , {-"t) *> 

(7) 



where the integral is over the BZ, the factor of 4 ac- 
counts for spin and valley degeneracies, and uj c (k z ) = 
eB/m*(k z ). After integrations over k p and the azimuthal 
angle at T = 0, the previous equation is reduced to 
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(8) 
where "^(O) is the zero-field conductivity, = k z c/2 and 



erB 
m*(0) 



(9) 



is the dimensionless parameter distinguishing between 
the regimes of classically weak (a <C 1) and strong 
(a 3> 1) magnetic fields. The limits of integration in 
Eq. pi are different for the electron and hole bands: for 
electrons, # m ; n = and max = 0q: f or holes , min = O 
and max = 7r/2, where O — cos~ 1 y> / eF/'2"/2 ~ 0.61 cor- 
responds to k z at which the Fermi energy intersects the 
bands (at k p = 0, cf. Fig. fl]). Likewise, the off-diagonal 
components are given by 
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The origin of the non-analytic dependence of the mag- 
netoconductivity on the magnetic field is evident already 
from Eqs. ffify and (10). Indeed, if the effective mass 
were independent of fc z , the weak-field behavior of a xx 
and a xy could be obtained by expanding Eqs. 



and 

( 10 1 in a. In our case, however, the region of integration 
over 9 for holes includes the H(H')- point (0 — tt/2), 
where the effective mass vanishes. An attempt to ex- 



pand Eqs. (|8|) and (|10|) in a leads to the (tt/2 - 0)- 1 
and ln(7r/2 — 0) divergences in a~ x and a~ y , correspond- 
ingly. Cutting off these divergences at n/2 — ~ a, 
we find that Aa xx (B) = a xx (B) — cr xx (0) oc \B\ and 
a xy (B) oc Sln|i3|. For electrons, the corresponding 
quantities are analytic: Aa+ X (B) oc B 2 and a xy (B) oc B. 
In the strong-field regime, a xx oc 1/B 2 and af y oc 1/B, 
as is expected from the Drude model. 



Integrals in Eqs. 
trary a. After some algebra 
the electron and hole contributions a^ = erf- 



and (10 1 can be solved for arbi- 
we obtain for the sum of 
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where a xx (0) = a xx + cr xx (0) = e ^f and j3 = 8/ [ir — 4# + 6sin(2#o)] = 1-26. For a <C 1, the equations above 
reduce to 



o xx {B) = a xx (0) [l - p {||a| - ? ya 2 + G(|a| 3 )} 
o xy {B) = -a xx (0)/3a 
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where 77 = 2 (tan 6> + ^ - 36» ) « 1.92. The lead- 
ing terms of these expansions coincide with the esti- 
mates obtained by cutting off the integrals in Eqs. ^ 
and ( |10p , and correspond to non-analytic field depen- 
dences of the magnetoresistivity: Ap xx (B) oc \B\ and 
p xy (B) oc -Bin \B\. Notice that p xy > for a«l, which 
indicates that weak-field magnetotransport is dominated 
by holes. Although we consider a perfectly compensated 
case, equal number densities of electrons and holes do 
not necessarily imply that a xy is equal to zero. This 
is so because the weak-filed limit of the Hall conduc- 
tivity in anisotropic conductors is not directly related 
to the volumes of the electron and hole Fermi surfaces 



but is given by an average of certain quantity (equal to 
v 2 /m*(k z ) in our case) over the Fermi surface. Although 
our system is compensated, the hole contribution to a xy 
exceeds that of electrons. In the strong-field limit, how- 
ever, the Hall conductivity must achieve a universal limit 
Cxy = ( n e ~ fl/O/e-B™ which is equal to zero for com- 
pensated metals (n e — nj/). Equation (lib I shows that 
this indeed the case: for a ^$> 1, the 1/B term in a xy (B) 
is absent and a xy (B) w — 2a xx (0)/3/15a 3 oc 1/B 3 , which 
is the expected behavior for compensated semimetals in 
the case when the 1/B 2 term is not allowed by lattice 
symmetry. At the same time, a xx (B) shows the usual 
behavior in the strong-field limit: a xx (B) oc 1/B 2 . 





FIG. 2. (Color online) Calculated dependence of 

Apxx/pxx{0) = [pxx (B) - pxx (0)] jpxx (0) on the magnetic 
field in graphite in the weak-f ield regime. Solid (blue): 
Apxx{B) corresponding to Eqs. |TTa| and ( |IIb[ ) for the non- 
relativistic spectrum, Eq. pi). Dotted (redj: Numerical re- 
sults for Ap X x(B) for the relativistic spectrum, Eq. p|. The 
dashed (black) line shows the asymptotic linear dependence 
at small fields: Ap xx /p xx (0) = Pna/2. Here a — er_B/m*(0). 



FIG. 3. (Color online) Calculated dependence of p xy /pxx(0) 
on the magnetic field in graphite in the weak-field regime. 
Solid (blue): p xy (B) corresponding to Eqs. (fTTajl and |TTb| for 
the non-relativistic spectrum, Eq. pi). Dotted (red): Numer- 
ical results for p xy (B) for the relativistic spectrum, Eq. p|. 
The dashed (black) line shows the asymptotic dependence as 

a -> 0: Pxy/Pxx{0) = f3a (in ^ - |V Here a = erB/m*(0). 



In passing we note that, although nonanalytic MR 
in graphite arises due to Dirac-like fermions near the 
H(H')- points, the effect occurs not because the spec- 
trum of these Dirac-like fermions is linear in k p , but be- 
cause the mass of normal massive fermions vanishes near 
the H(H')- points rendering the fermions Dirac-like. In 
fact, o~ xx of strictly two-dimensional Dirac fermions in 
doped graphene is expected to behave as B 2 for T <C Ef- 
(A finite-temperature correction to a xx and the leading 
term in p xx scale as \/B but are both exponentially small 
in this temperature regime!^) 



B. "Relativistic" effects in the energy spectrum: 
a new scale for magnetic field 

In the previous Section, we found that the magne- 
toconductivity behaves nonanalytically as a function of 
the magnetic field in the weak-field regime, when trans- 
port is controlled by extremely light holes residing near 
the H (_ff')-points. This result, however, is valid only in 
the non-relativistic approximation for the energy spec- 
trum, Eq. (KM), which is obtained by expanding the rel- 
ativistic spectrum, Eq. pi) in k p . Returning to Eq. pi), 
we see that the actual expansion parameter is the ra- 
tio 2v p kp/"/iT. The width of the regions contributing 
to non-analyticities in the magnetoconductivty, Sk z — 
tt/c — k z ~ a/c, shrinks in proportion to the magnetic 
field. Since T « c5k z /2 ~ a < 1 near the H(H')- 
points, we must eventually reach such a weak magnetic 
field when 2v p k p /jiT ~ 1, and thus the expansion breaks 



down. Using typical k p estimated as \ef\ 



'v p for Fermi- 



surface carriers, we find that the expansion is valid only 
for a a = |e_f|/7i <a<l. In very weak magnetic fields, 
such that a < |ef|/7i ~ 0.06, one needs to consider the 



complete spectrum in Eq. pi) and recalculate the depen- 
dence of the conductivity on the magnetic field. It is 
straightforward to show that in this limit the magneto- 
conductivity is again analytic, i.e., Aa xx (B) ex B 2 and 
o~ xy (B) ex B. Since spectrum pi is isotropic in the ab 
plane, the Boltzmann equation can still be solved ana- 
lytically. However, the resulting integrals for the mag- 
netoconductivity need to be solved numerically. Figures 
[2] and [3] show comparisons of the magneto- and Hall re- 
sistivities in the relativistic and non-relativistic models. 
At weak fields, the curves for both the relativistic as well 
as the non-relativistic models show similar non-analytic 
dependences. For example, fitting Ap xx / p xx (0) into a 
linear function of B in the weak-field region (we take 
0.06 < a < 0.25 since the linearity is most pronounced 
here) yields slopes that differ by w 30%. However, the 
absolute values of Ap xx /p xx (0) in the two models are 
different: e.g., at a = 0.2, the two results differ by 
f=a 60%. This difference in the values is due to the fact 
that at superweak fields the linear field dependence of 
the magnetoresistivity predicted by the non-relativistic 
model is replaced by an analytic, quadratic field depen- 
dence predicted by the relaltivistic model thus creating 
an offset. Note that the mechanism described above also 
regularizes a logarithmic divergence of the Hall constant 
Rh = p xy /B ex ln|B| implicit in Eq. (lib). 



To summarize, we see that, contrary to the case of con- 
ventional metals, the magnetoconductivity of graphite 
exhibits three rather than two characteristic regimes: 
superweak, weak, and strong magnetic fields. In the 
superweak-field regime (a ^C ao), the magnetoconduc- 
tivity is analytic: Aa xx (B) ex B 2 and a xy ex B. In the 
weak- field regime (ao <Ca€ 1), both components of a 
are non-analytic: Aa xx (B) ex \B\ while a xy ex _Bln|i?|. 
In the strong-field regime (a 3> 1), the magnetoconduc- 



tivity behaves as is expected for a compensated semi- 
metal. The occurrence of a new scale for the magnetic 
field, ao> is due to a peculiar feature of the graphite en- 
ergy spectrum: the existence of Dirac-like fcrmions near 
the H(H')- points. Notice that linear MR, which is the 
central result of this section, is still a classically weak-field 
phenomenon. This should be contrasted to other pro- 
posed mechanisms (either classical or quantum), where 
linear MR occurs in strong fields.^ 



Discussion 
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1. Macroscopic inhomogeneities 

In the previous section, we have shown that the spec- 
trum of graphite allows for linear MR in a certain interval 
of the magnetic field. Using typical values for r and band 
parameters for graphitepESJ we find that the range of lin- 
ear MR, ao < a < 1, translates into a range of fields 
from ~ 0.006 T to ~ 0.1 T. In experimentjinear MR is 
indeed observed in some graphite samples Q2HHI however, 
it is found to exist not only in this field range but in 
stronger fields too, which cannot be accounted for by a 
simple model of graphite employed in the previous sec- 
tion. Therefore, other factors are also probably at play. 

Linear MR in classically strong fields is often ascribed 
to macroscopic inhomogeneities in the sample. If the 
charge carriers in the material, instead of being uniformly 
distributed, form macroscopic puddles with different Hall 
conductivities, then the effective magnetoresistivity of 
the sample is linear in classically strong fields. This 
idea has been explored in the past pS and recent observa- 
tions of linear MR in some other materials, such as silver 
chalcogenides, have been attributed to the presence of 
such inhomogeneities J3SEU It is plausible that linear MR 
in graphite in strong fields also originates from such inho- 
mogeneities. A general case of unequal volume fract ions 
has been studied under various approximations^!^ hut, 
just to illustrate the point, we use the exact resuHpSfor 
the case of a two-dimensional material with equal par- 
tial volumes. (Strictly speaking, the material can still 
be three-dimensional, only the inhomogeneity has to be 
two-dimensional) . The effective conductivity a e of such 
a system can be written as 



a ( V 2) 

u XX v XX 




1/2 



a (2) a (1) 

u xy v xx 



-CT ( V 2) 
u xy u xx 



U XX H^ U XX 



(13a) 
(13b) 



where a^ (i = 1,2) are the conductivity tensors of the 
individual components. 
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It is easy to see from Eqs. (|13a| and (|13b) how MR 

1/B 2 and 



x 



becomes linear in strong fields: 

Oxy oc 1/B in this regime, then a xx oc 1/B instead of be- 



FIG. 4. (Color online) Calculated dependence of 

Ap xx /p xx (0) on the magnetic field in graphite over a wide 
range of fields. Solid (red): an inhomogeneous system de- 
scribed by Eqs. (13a l and (13b) . Dashed (blue) and dotted 
(black): compensated and decompensated homogeneous sys- 
tems, respectively. Here a = erB/m* (0) . 



ing proportional to 1/B 2 as in the homogeneous case. On 
the other hand, since a x behaves in the usual way, i.e., 

as 1/B, we have p xx = ^f} - ^ oc B in strong fields. 

,(0 



("•:.)'+(« 



For graphite, the required 1/B dependence of a X y implies 
that each of the components must be decompensated 
Decompensation can be modeled by shifting the Fermi 
energy away from its value for a compensated system 
We define the degree of decompensation as £ = Sef/sf 
and choose £ = 0.09 and £ = —0.01 for components 1 and 
2, correspondingly. Figure H shows Ap xx [solid (red)] for 
an inhomogeneous system. The magnetoconductivities 
of the individual components were calculated using the 
relativistic spectrum from Eq. (pi). Also shown is Ap xx 
for a compensated [dashed (blue)] and decompensated 
[dotted (black)] system. (We used £ = 0.05 for a homo- 
geneous decompensated system.) Comparing the solid 
and dashed curves, we see how inhomogeneity transforms 
quadratic MR into linear one. The dotted line illustrates 
that decompensation leads to saturation of MR in a ho- 
mogeneous system. 

It is thus possible that linear MR observed in strong 
fields arises due to extrinsic effects, such as macroscopic 
inhomogeneities. This effect is distinct from intrinsic lin- 
earity discussed in Sec. |IH A| which arises due to the 
presence of Dirac-like holes in graphite and is a weak- 
field phenomenon. Note that intrinsic linearity is not 
affected by decompensation: the effect of the latter be- 
comes pronounced only in strong fields where it leads to 
saturation of MR. 



2. Quantum magnetooscillations 

Quantum magnetooscillations arise on top of classical 
MR in strong magnetic fields when w c is larger than not 
only 1/r but also temperature. Magnetooscillations owe 



their origin to the extremal cross sections of the Fermi 
surface perpendicular to the field. In graphite, the ex- 
tremal (maximal) cross-sections of both the electron and 
hole Fermi surfaces occur away from the H(H')- points. 
The massive carriers located at these cross-sections give 
rise to two sets oscillations in the resistivity. However, at 
the H(H')- points where the BZ ends in the z-direction, 
the Fermi surface remains open, leading to tiny but non- 
zero (minimal) cross-sections (cf. Fig. fTl). Therefore, 
it is reasonable to ask if quantum oscillations can also 
arise due to these minimal cross-sections harboring Dirac 
fermionsPIS We answer this question in the negative: 
within the accepted band structure, it is not possible 
to have a new set of oscillation arising from the Dirac 
fermions. To see this, it suffices to recall that the oscilla- 
tory part of the conductivity leading to the Shubnikov-de 
Haas effect (or of magnetization in case of the de Haas- 
van Alphen effect) is given by 



q^O 



dk z I n (k z )e i2 ^ k »\ 



(14) 



where I n (k z ) is a combination of the single-particle 
Green's functions depending on the quantity being cal- 
culated and n{k z ) is the Landau index as a function 
of k z . In the limit of a large number of Landau levels 
(n 3> 1), one computes the integral via the stationary 
phase approximation,^ in which the main contribution 
to the result comes from those values of k z where n(k z ) 
has an extremum. In case of graphite, the Landau levels 
with index number n for the Dirac-like Fermions near the 
H(H')- points are described by^l 



1 



2ev 2 p B 



e n < £„ - A - 271 cos 



(15) 



It is obvious that this expression does not have a non-zero 
derivative, and hence an extremum value, at the H(H')- 
points where k z = ±tt/c. Therefore, no new quantum 
oscillations are expected to result from the carriers near 
the H(H')- points. Note that the sole reason for this 
negative result is an extra factor of 2 in the denominator 
within the cosine term in Eq. ( 15 ). This factor of 2 results 



from Bernal-stacking of the graphene layers in graphite 
which makes the periodicity in the z-direction to be two 
lattice planes instead of one. If one were to construct 
artificial graphite by placing graphene layers directly on 
top of each other in the AAAA... fashion, Dirac fermions 
at the H(H')- points would indeed give rise to their own 
set of quantum oscillations. Note also that our argument 
means only that Dirac fermions do not give rise to new 
oscillation frequencies. The issue of the oscillation phase 
is beyond the leading-order semiclassical approximation 
employed here. 



IV. CONCLUSIONS 

In conclusion, we have shown that there exists an in- 
terval of the magnetic fields in which in-plane magne- 
toresistance in graphite (with B along the c-axis and the 
current in the a6-plane) scales linearly with the field. Ex- 
tremely light, Dirac-like carriers located near the H(H')- 
points contribute to such a non-analytic behavior. Lin- 
ear magnetoresistance occurs for classicaly weak fields, 
unlike other mechanisms where such a behavior is found 
in classically strong or even ultraquantum fields. The 
Hall resistivity also show a non-analytic, B In \B\ field de- 
pendence and is dominated by holes. However, observed 
linear magnetoresistance sometimes spans over the en- 
tire range- from weak to classically strong fields, and 
beyond - which cannot be accounted for by the mecha- 
nism described above and could be due to macroscopic 
inhomogeneities in real samples. The light carriers near 
the H(H')- points, however, do not give rise to quan- 
tum oscillations of their own, the reason being the Bernal 
stacking of graphene planes in graphite. 
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